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Abstract. We investigate the caustic structure of a lens 
composed by a discrete number of point -masses, having 
mutual distances smaller than the Einstein radius of the 
total mass of the system. Along with the main critical 
curve, it is known that the lens map is characterized by 
secondary critical curves producing small caustics far from 
the lens system. By exploiting perturbative methods, we 
derive the number, the position, the shape, the cusps and 
the area of these caustics for an arbitrary number of close 
multiple lenses. Very interesting geometries are created in 
some particular cases. Finally we review the binary lens 
case where our formulae assume a simple form. 

Key words: Gravitational lensing; Close binary stars; 
Clusters of galaxies; Quasars 



1. Introduction 

The binary Schwarzschild lens is one of the most inten- 
sively studied model. In fact, in a relatively simple way, 
it shows many features that are observed in general grav- 
itational lenses, such as the formation of multiple images, 
giant arcs and a not trivial critical behaviour. The first 
study about the binary lens with equal masses was made 
by Schneider & Weifi (1986). They derived the critical 
curves and the caustics showing that three possible topolo- 
gies are present depending on the distance between the 
two lenses. Erdl & Schneider (1993) extended these results 
to a generic mass ratio of the two lenses. Witt & Petters 
(1993) reached the same results using complex notation. 
In some limits, Dominik (1999) enlightened the connection 
between the caustics of the binary lens and other models, 
such as the Chang-Refsdal lens (Chang & Refsdal 1979; 
1984) and the quadrupole lens. 

The critical curves and the caustics of multiple lenses 
can develop very complicated structures, so that the at- 
tempts to gain some information about them have been 
very few. However there is a great interest in this prob- 
lem for its applications in particular situations, such as 
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planetary systems (Gaudi et al. 1998), rich clusters of 
galaxies and microlensing of quasars by individual stars 
in the haloes of the lensing galaxies (Chang & Refsdal 
1979; Kayser et al. 1988). 

In some special situations, the critical curves of multi- 
ple lenses can be derived by perturbative methods, refer- 
ring to the single Schwarzschild lens as the starting point 
for series expansions (Bozza 1999; Bozza 2000). These 
methods work very well in planetary systems, for a lens 
very far from the others and systems where mutual dis- 
tances are very small with respect to the total Einstein ra- 
dius. In the first two cases, the complete caustic structure 
has been derived and the connections with other mod- 
els have been showed. In the last case, only the central 
caustic coming up from the deformation of the total Ein- 
stein ring has been studied. Besides this main curve, there 
are many small critical curves forming among the masses. 
The caustics generated by these curves generally lie far 
from the centre of mass and can have some influence on 
sources distant from the mass distribution. Moreover, they 
move very quickly as the parameters of the system change 
(Schramm et al. 1993) constituting the most problematic 
feature to control in numerical simulations. For these rea- 
sons they are sometimes dubbed ghost caustics. In rapidly 
rotating binaries they may have superluminal projected 
motion requiring a non-static treatment of light deflec- 
tion (Zheng & Gould 2000). 

In this paper we use complex notation to face the prob- 
lem of secondary caustics of close multiple lenses. In this 
way we can study them as deeply as the other caustics, 
completing the previous works. We shall see that different 
classes of secondary caustics can be recognized, showing 
different geometries. 

After some review of multiple lensing in Sect. 2, in 
Sect. 3 we calculate the number and the position of sec- 
ondary critical curves for an arbitrary number and con- 
figuration of lenses. Then, in Sect. 4, we treat the simple 
caustics and in Sect. 5 the multiple caustics (the distinc- 
tion will be explained at the end of Sect. 3). In Sect. 6 we 
specify our formulae for the binary case and in Sect. 7 we 
give the summary. 
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2. Basics of multiple lensing 

We shall study a system of n point -lenses placed at po- 
sitions Xi = {xi\]Xi2) in coordinates normalized to the 
Einstein radius 



R-v. 
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OL 



Dos 



(1) 



where M is a reference mass (it can be chosen to be the 
total mass, the typical mass of a single object or anything 
else). The source coordinates y = (yi;y2) are normalized 
to the scaled Einstein radius -Rft-R 2 ^- . The masses m, of 
the lenses are measured in terms of Mq. 

We introduce the complex coordinate in the lens plane 
z = x\ + ix2 and the complex source coordinate y = ij\ + 
iyi- The positions of the masses will be denoted by z% = 
xn + ixi2- We also introduce the functions 



~t - z i) 



(2) 



The lens equation for our system of n masses reads 
(Witt 1990) 



Si (z) 



(3) 



Given a source at position y, the z's solving this equation 
are the images produced by gravitational lensing. 

This map is locally invertible where the determinant 
of the Jacobian matrix 



det J =1 — 



dydy 
&z &z 



= i-is 2 (*)r 



(4) 



is different from zero. The points where the Jacobian de- 
terminant vanishes are arranged in smooth closed curves 
called critical curves. The images of these points through 
the lens map (j^) in the source plane are called caustics. 
When a source crosses a caustic, creation or destruction 
of pairs of images occurs and the magnification diverges 
(Schneider, Ehlers & Falco 1992). 

This is all we need to start our search for secondary 
caustics in close multiple systems. The fundamental hy- 
pothesis we make is 



(5) 



\zi\ < VM Vi 



where M is the total mass of the system. In this way, the 
distances between pairs of lenses will be very small with 
respect to the Einstein radius of the lens that we would 
have if all the masses were concentrated at the origin. This 
Einstein radius is \[M in our notation. The relation (||) 
allows us to consider the S clS pcrturbative parameters in 
a series expansion. Then we can solve the equation det J = 
at each order, writing its solutions as series expansions 
in powers of the perturbative parameters. In this way we 
shall find the critical curves of this system and study their 
properties analytically. 



3. Number and positions of secondary critical 
curves 

Close multiple lenses have two classes of critical curves: 
the main critical curve, resulting from the deformation 
of the Einstein ring of the total mass lens, and the sec- 
ondary critical curves, forming inside the distribution of 
the masses. 

Effectively, if we multiply the equation det J = by 



the quantity Yi \ z ~ z i 

i=i 

and ~2\ 



we get a complex equation in z 



m 



= 0. 



(6) 



At the zero order, putting all Zj's to zero, this equation 
becomes 



\zf n -*(\z\ 4 



M 2 = 0. 



(7) 



This equation has the solution \z\ = VM, that is the Ein- 
stein ring of the total mass lens. Taking this solution as 
the starting point of a perturbative expansion, we get the 
main caustic. The details of this calculation are in (Bozza 
2000). But the presence of the solution z = indicates 
that also this value can be taken as the starting point for 
another expansion. This is just the value we shall take to 
find the secondary critical curves. 

Having observed the zero order situation, we can start 
our perturbative approach, searching for the first order so- 
lution. Then we write the solution z as a series expansion: 



Z = Z +o(\Zi\) , 



(8) 



where zq is of the first order in \z%\. Stopping at the first 
order, we put z = zq in Eq. (^j). We see that the first 
term becomes of order 4n, while the second is of order 
An — 4. Then the latter dominates the first and Eq. (^|) is 
equivalent to 



i=l 



3#« 



0. 



(9) 



This is a polynomial equation of degree 2n — 2. Then, 
for a system of n close lenses, there are, at most, 2n — 2 
points where the Jacobian determinant vanishes (at the 
first order in Zi), corresponding to 2n — 2 secondary crit- 
ical curves. This is the first main result of our work. It is 
consistent with the binary lens, since two secondary criti- 
cal curves are predicted by this formula. 

Eq. (^|) can be solved analytically for two and three 
lenses, otherwise we have to resort to simple numerical 
methods. In Sect. 6, we shall specify these and the fol- 
lowing results for the binary lens where a manageable ex- 
pression for the positions of the critical curves is available. 
For the triple lens, the analytical solutions are too cum- 
bersome to allow a detailed study. 
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Now, we have a straightforward way to calculate the 
positions of the secondary critical curves for an arbitrary 
configuration of close multiple lenses. Then, we can avoid 
the traditional blind sampling of the Jacobian determinant 
on the lens plane and reach, by this new method, the full 
efficiency. 

We take the generical solution zq of Eq. (^) as the first 
order term of our expansion. From now on, we use the 
notation 



The first row is the order —2 and is null according to Eq. 
(pi]). The second row is the order —1 and the third row is 
the order zero. Inserting this expansion in (|l4j), the lowest 
order equation is of order —2: 



M 2 |s 3 °| 



0, 



(16) 



Being zq a simple root, S3 ^ 0, so that £2 = 0. 

The successive terms in the expansion of the equation 
(O) are of order zero: 



s k = S k ( z o) 



(10) 



1 



I 9° I 



0. 



As both z and Zi are of the first order, according to our From this equation we have 
perturbative expansion, S° has all denominators of order 



k and then it is of order — k. 

To continue our study we do not need an analytical 
expression for z§. We shall just use the fact that zq is a 
solution of Eq. (Q), that is equivalent to say that 



S 2 ° = 0. 



(11) 



Of course, we have to distinguish between simple roots 
of Eq. (^) and roots of higher multiplicity. Remembering 
that the k th derivative of S 2 (z) is proportional to Sk+2 (z), 
we have the equivalence between the following statements: 

zq is a root of multiplicity p <^> S° +2 =0 Vfc < p. (12) 

We shall treat separately the caustics coming from sim- 
ple roots (hereafter called simple caustics) and the caustics 
coming from multiple roots (hereafter multiple caustics). 



4. Simple caustics 

These caustics are largely the most common as we explain 
in the next section. So they surely have the most practical 
interest. 



4.1. Shape of the critical curves 

Once found the positions of the critical curves, we can 
carry on our perturbative expansion to discover the shape 
of these curves. So we put 



zq + (-2 + £3 



(13) 



where zq is the position of one simple critical curve, found 
by Eq. (^|), and ej are the corrections of order \zi\ 3 . It is 
convenient to use the original equation det J = 0, which 
can be written in the form 

1 - S 2 (z + e 2 + e 3 ) S 2 [z + e 2 + e 3 ) = 0, (14) 

starting from Eq. (^) 

The expansion of S 2 is 

= S°2 + 

-2e 2 S 3 ° + 

-2e 3 Si + 3ejSl + ... . (15) 



S 2 {zq + e 2 + e 3 ) 



2|S 3 °| 



(17) 



(18) 



Then the third order contains the first information on the 
shape of the critical curve. Eq. ( fl8l ) tells us that the critical 
curve at position zq is a circle centered on zq with radius 

r =mi (19) 

From the form of S3 (see Eq. (|J)), we see that the closer 
the critical curve is to some mass, the higher the value of 
S3 , the smaller the radius of the circle. 

If we multiply all masses by a factor A, the positions 
of the critical curves do not change, because A factors 
out from Eq. (^|), but their radii change as A -1 . If we do 
the same with the positions of the masses instead, the 
positions of the critical curves scale as A and their radii 
scale as A 3 . 

4-2. Caustics 

To find the caustics corresponding to the simple critical 
curves, we just have to put the critical curve, in its obvious 
parameterization 

z (9) =z a + re 10 < 9 < 2tt, (20) 

into the lens equation (^|) and expand to the third order: 



-2i6 



s° 3 



(21) 



We can observe that the lowest order is —1 and is in- 
dependent on 9. It represents the position of the caustic. 
From the order of this term, we can deduce that these 
caustics can lie very far from the origin of our system, 
going to infinity as the distances among the masses are 
reduced to zero. The successive term is z , which is of the 
first order and represents a correction to the position. Fi- 
nally, the shape of the caustic is given by the third order. 

The cusps of a caustic are characterized by the vanish- 
ing of the tangent vector. To find them, we have to require 
that 



dy(9) 
d9 



= 



(22) 
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and solve for 9. Taking y (9) from Eq. (^l|), this equation 
can be simplified into 




whose solutions are 
0k = -\* 



^ (s§) 



2kn 



0,1,2, 



(23) 



(24) 



where arg yields the argument of a complex number. 

We have three cusps. So, in any close multiple system, 
having only simple secondary caustics, these caustics have 
a triangular shape. 

Finally, we calculate the area of these caustics. This 
can be done by the integral 



.4 



2/2dyi 



(25) 



where 7 is the caustic in its clockwise direction. We have 
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A = ~J[y{9)-y (6)] d e [y (9) + y (9)] d9. 



(26) 



The minus in the right member comes from the fact 
that our parameterization is counterclockwise. The inte- 
gral only involves complex exponential functions and the 
result is 

A=\irr\ (27) 

So the extension of the simple caustics is of the sixth order 
in the separations among the lenses, justifying the evasive 
nature of these caustics. 

With our expansions, we have attained considerable 
analytical information about the secondary caustics es- 
tablishing their shape, the area, the number of cusps in a 
completely general way. However, since these results are 
the fruit of perturbative approximations, it is important 
to discuss their accuracy. So we propose a comparison be- 
tween our perturbative results and the numerical ones in 
a typical situation. We consider a system constituted by 
three lenses disposed as in Fig. |l|. According to our previ- 
ous statement, this system can form, at most, four simple 
secondary critical curves. For our choice of parameters, 
we display their positions in the same figure. The caus- 
tics produced by these curves are shown in Fig. 2 where 
they are compared to the numerical ones. We have taken 
the distances among the masses of this distribution to be 
one tenth of the total Einstein radius. Even for this not 
too small value, the positions and the shapes of the sec- 
ondary caustics are reproduced with a striking accuracy. 
It is also to be noted that the quality of numerical results 
is improved thanks to the guide provided by perturbative 
results. 

So we see that the analytical formulae derived in this 
section are very good approximations to the quantitative 
characteristics of the secondary caustics, proving to be 
highly reliable. 



Fig. 1. Positions of the four secondary critical curves (in- 
dicated by the small dots) for a triple system with masses 
mi = 0.25, m2 = 0.25, 7713 = 0.5 and positions z\ = 0.1, 
z 2 = -0.1, z 3 = 0.014 + i0.089 (indicated by the three 
crosses). 



5. Multiple caustics 

In this section, we consider the case where zq is a mul- 
tiple root of Eq. (^|). The parameters space of a system 
with n lenses is 3n — 4 dimensional, since each mass adds 
three parameters (its mass and its coordinates in the lens 
plane). Four parameters can be eliminated by considering 
equivalent those systems differing by a global translation 
or rotation and/or by a global scale factor. Thus, for ex- 
ample, the binary lens is completely characterized by the 
mass ratio and the separation between the lenses. 

The requirement of a double root in Eq. (|^) translates 
into the vanishing of the derivative of this equation with 
respect to z. This is one constraint equation, then the 
points of the parameters space producing multiple roots 
constitute a 3n — 5 dimensional hypersurface, thus having 
measure zero. For this reason, the occurrence of multi- 
ple roots is relatively rare. Anyway, very interesting fea- 
tures emerge, justifying a detailed study of these particu- 
lar cases. 



5.1. Critical curves 

Suppose that zq is a root with multiplicity p. We have to 
find the correct order of the perturbation to insert in the 
equation det J = 0, representing the shape of our critical 
curve. According to the equivalence (|l^), the S , " +2 with 
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Fig. 2. Here are the four caustics 
produced by the critical curves of 
Fig. [j] from left to right. The solid 
curve is the perturbative caustic and 
the dashed line is the numerical one. 



k < p are null. Then, we put 

z = z + e, 



(28) 



with the order (that we shall indicate by q) of e to be 
found. We only assume that q be higher than one. Then, 
the expansion of S 2 {%) is 



S 2 (z + e) = 5° - 2eS$ + SeS^ 



+ (-ir(k + l)e k S° k+2 + 



(29) 



The fc th term is of order qk — (k + 2) and the first 
term to be non-null is that for k — p. When we put this 
expansion in the equation det J — 0, the first non-null 
term is 



-{p+l) 2 \^S p+2 \ 2 



(30) 



having order 2qp — 2 (p + 2). If this order is less than zero, 
we just get from det J = that e = 0, but if the order of 
this term is zero, then the zero order expansion of det J = 
also involves another term (equal to 1): 



l-(p+l) 2 \ePS p+2 \ 2 = 
and the equation gives the non-trivial solution 

1 



[(P+1)I3 



p+2\ 



(31) 



(32) 



This happens when the order of e is q — 2±2. This is con- 
sistent with the result of the previous section, because, 
for p = 1, q = 3. For p = 2, we have that the first non 
trivial order is the second and, for p = 3, it is the order 
|. When p increases, the order of this perturbation de- 
creases, approaching 1 as a limit. This means that at high 
multiplicities, the perturbative expansion becomes always 
less accurate, requiring ever more terms for an adequate 
description of the caustics. Anyway, the main characteris- 
tics of the caustics can be derived retaining just the first 
correction and that is what we shall do. 

The critical curve just derived is again a circle with 
radius 



1 



[( P +l)\S p+2 \] 1/p ' 
becoming greater with the multiplicity. 

5.2. Caustics 

We take, as before, the parameterization 
z (6) =z Q + re 19 



(33) 



(34) 



for the critical curve, with r given by Eq. (|33|). Putting 
this expression into the lens equation (|^) and expanding 



Valerio Bozza: Secondary caustics in close multiple lenses 



o 




Fig. 3. Critical curves in the event of a double root. The 
masses are indicated by the crosses. The double critical 
curve is the one at the bottom-center, while the other two 
are on the top-left and top-right and are hardly visible in 
this picture. 



to the q order, we get 



y (9) = -S 1 + z + r 



e w + (-If ■ 



7 p+2 



J p+2 



(35) 



The q th order is the first depending on 9 and deter- 
mines the shape of the caustic. 

To understand this shape, we calculate the cusps as in 
the previous section. The equation for the cusps is 



e (p+2)*e 

and its solutions are 
a (-1) P 



+ (-1) 



P+i 



°p+2 



= 



(36) 



p+2 



P 



g (s° P - 



P+2) 



2kn 
p + 2 



0<k<p+2. (37) 



Now we have p + 2 cusps. This is a very interesting 
result, because the caustic assumes the shape of a regular 
polygon with p + 2 curved sides. 

The area of the multiple caustic can be calculated in 
the same way as for the simple one. We just give the result: 



A = 



P 2 

-irr . 



(38) 



Fig. 4. Caustic corresponding to the double critical curve 
of Fig. ||. The dashed curve is the numerical caustic and 
the solid line is the perturbative one. 



It is of order 2^±£. So, increasing the multiplicity from 1 
to infinity the order of the area lowers from 6 to 2 and the 
extension of the caustic becomes ever more important. 

Some other consideration about the limit for p — ► oo 
can be done. The number of cusps become infinite and, 
from Eq. (|35|), we see that the caustic becomes a circle of 
radius r. In fact, the area becomes ttt 2 . 

5.3. An example: a double caustic in a triple lens 

Now we shall practically see how our formulae work in 
the case of a multiple caustic. We consider three masses: 
mi = 0.25, m2 = 0.25 and 7773 = 0.5. We fix the positions 
of the first two: z\ = 0.1, Z2 = —0.1; but we let the third 
free for the moment. We simultaneously solve Eq. (||) and 
its derivative with respect to zq, for the two unknowns zq 
and 03. We find six possible positions of the third mass, 
giving rise to a double root of the positions equation. None 
of them is a triple root. Two of these positions are on 
the £2-axis. We choose one of them: z% = iO. 084263. The 
double root is in z = — 7,0.0299. In Fig. ||, we see that 
the critical curve in this point is much greater than the 
other two. In fact, the radius of the double critical curve, 
calculated by Eq. (33), is 8.49 x 10~ 3 , while the radius of 
the two simple critical curves is 6.13 x 10~ 4 , according to 
Eq. ©. 



In Fig. |4|, we show the caustic generated by this double 
critical curve. The geometry is correctly predicted by our 
perturbative expansion: there are four cusps in a double 
caustic. We see that the approximation is less accurate 
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than before, as we anticipated in our discussion about the 
order of the perturbation. However, for double caustics, it 
is not so difficult to add another term to the perturbative 
expansion and reach the same accuracy of the simple caus- 
tics. The third order term in the critical curve depends on 



e 3 = 6r 2 Rc 



The successive term in the caustic is 

-3i9 2 



e™e 3 + 3e~ MO r z e 3 sl 



r 4 e- 4 * e sl 



(39) 



(40) 



Double caustics, and, more generally, multiple caus- 
tics, are formed by the union of small caustics, in some 
sense. Another interesting question is: what happens if we 
change the parameters in the neighbourhood of our par- 
ticular choice producing the double root? We expect the 
double critical curve to separate into two smaller ovals 
and the double quadrangular caustic to break into two 
triangular ones; but this can happen in different ways. 

In this regime, the perturbative caustics are simple. 
However, as the parameters tend to give the double root, 
S3 tends to zero, yielding a diverging r for the simple crit- 
ical curves, according to Eq. ( [l9| ) . The transition with the 
formation of the double critical curve is thus not repro- 
duced. Guided by perturbative approximations, the break 
of the double caustic, when z 3 moves out from the position 
iO. 084263, can be investigated numerically. The results are 
shown in Fig. |^. 

In case a, z 3 = iO. 082787, i.e. we have moved the third 
mass towards the others. The critical curve breaks in the 
horizontal direction. Looking just at the thick line in Fig. 
||a2, representing the numerical caustic, we see that the 
top cusp and the bottom cusp develop a butterfly geome- 
try. At some critical value, these butterflies touch and the 
two resulting triangular caustics move away along the hor- 
izontal direction. We have displayed in the same plot the 
perturbative caustics too. Obviously, they are simple caus- 
tics, so they cannot show the butterfly geometry but they 
can help in understanding how the separation occurs. We 
also notice that the simple caustics cover the area of the 
numerical transition double caustic very well constituting 
a significant approximation anyway. 

In case b, z 3 = iO. 085759, so that the third mass is 
farther from the others. Now the critical curve breaks in 
the vertical direction and so does the caustic. The left and 
the right cusps transform into butterflies. These butterflies 
are slightly distorted by the fact that the resulting simple 
caustics have different sizes: the one on the top is smaller 
than the other. 

In case c, z 3 = 0.00147 + ^0.084263. We have displaced 
the third mass in the horizontal direction. The critical 
curve breaks diagonally and so does the caustic. But this 
time the transition occurs with a simple beak-to-beak sin- 
gularity rather than with butterflies. While in the previous 
situations the two simple caustics in the last step of the 
separation touch with a fold, here they touch with a cusp. 



6. Secondary caustics in binary lensing 

In this section, we specify our results for the binary case 
where simple analytical formulae can be written. Let's 
consider two masses placed on the horizontal axis and let's 
choose the origin in the centre of mass. We call the separa- 
tion between the masses a, then we have z\ — m f a and 
mia . Eq. <M) is of second degree. Its solutions 



Z2 - 

are 



mi+tii2 ' 



Z0 



rri2 — mi ± iJmim2 

1 . 

mi + TO2 



(41) 



They are always simple and lie on a circle of radius a/2 
centered in the middle of the two masses. 
The radius of the two ovals is the same: 



^Jmim2a 
2 (mi 



m 2 ) 



(42) 



Its maximum value 



4M„ 



is reached when the two masses 



are equal, in fact, in this case, their distance from the two 
masses is maximum. 

The two caustics are given by the following expression: 



m\ -mi 2Jmim 2 
V (<?) = + 1— 



m2 — 77ii , . Jmim 2 
a ± la 

Till + 7712 7711 + 7712 



+ a ^/ m i?ii2 



,.i0 



± 



te 



-2i8 



2(7711+7712) A(y/m^±i v ffn2~) 
Their cusps are at positions 



(43) 



arg 



Their area is 



2kn 



k = 0,1,2. 



7TT71l77l2a" 

; (7711 + m 2 ) 4 



(44) 



(45) 



reaching the maximum value . m the equal masses 
case. 

7. Summary 

Multiple Schwarzschild lensing represents a very rich ter- 
rain for the exploration of caustics in gravitational lensing. 
The occurrence of different kinds of singularities stimu- 
lates new investigations. 

In this paper we have applied perturbative methods 
to secondary caustics, forming when the masses are close 
each other with respect to the total Einstein radius. In 
this way we have been able to establish the number of the 
caustics for any lens configuration, the positions and the 
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Fig. 5. Critical curves and caustics for z 3 close to iO. 084263. The left column shows the critical 
curves and the right column the caustics. The thick lines are the numerical curves and the thin 
lines are the perturbative ones. The choice of zz in cases a, b and c are given in the text. 
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shapes, with a complete characterization of the geome- 
tries arising in all cases. Moreover, quantitative formulae 
for the area and other features of these objects have been 
given. As we have seen, in the most common case, the 
shape of the simple caustics is always triangular. Anyway, 
multiple caustics exist developing a great variety of be- 
haviours, giving rise, to curves having a number of cusps 
ranging from four to infinity. The breaking of multiple 
caustics can follow different ways depending on how the 
parameters of the system change. 
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